Low-temperature properties of Pb(Zri_a;Ti^)03 solid solutions 
near the morphotropic phase boundary 
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A first-principles-derived approach is used to study structural, piezoelectric and dielectric properties 
of Pb(Zri_a;Ti2,)03 (PZT) solid solutions near the morphotropic phase boundary at low tempera- 
ture. Three ferroelectric phases are found to exist: a tetragonal phase for larger x compositions, 
a rhombohedral phase for smaller x compositions, and the recently discovered monoclinic phase in 
between. In this monoclinic phase, the polarization associated with the Zr atoms behaves differently 
from the polarization associated with the Ti atoms. As the composition x decreases, the former 
rotates more quickly towards the pseudo-cubic [111] direction and grows in magnitude, while the 
latter lags in its rotation and its magnitude shrinks. The local microscopic structure is found to 
deviate significantly from the average structure in these PZT alloy phases as a result of fluctuations 
in the directions and magnitudes of the local polarizations. The monoclinic phase is characterized 
by a very large piezoelectric and dielectric response. 
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I. INTRODUCTION 

Since the beginning of the 70's, Pb(Zri_a;Tia;)03 
(PZT) alloys have been known to exhibit a morphotropic 
phase boundary (MPB) separating a ferroelectric region 
with a tetragonal ground state {x > 0.52) from a ferr 
electric region with rhombohedral symmetry (x < 0.45)1 
Very recently, synchrotron x-ray powder diffraction stud- 
ies have revealed that there is in fact a third phase in the 
vicinity of the MPB at low temperaturea. This phase is 
ferroelectric, adopts a monoclinic symmetry and occurs 
within a narrow range of composition x that separat' 
the tetragonal region from the rhombohedral regionB' 
It has been suggested that this monoclinic phase acts as 
a transitional bridge between the tetragonal phase, for 
which the electrical polarization P lies along the pseudo- 
cubic [001] direction, and the rhombohedral phase, for 
which P is along the pseudo-cubic [111] directioncl. The 
use of a recently developed first-principles approach has 
confirmed that this is indeed the case, since the polar- 
ization of the monoclinic phase has been found to rotate 
from the [001] to the [111] direction as the composition 
X decreases in the MPB regions. 

Despite these experimental and theoretical advances, 
some features related to this MPB are still unknown. For 
instance, one may wonder what is the separate contribu- 
tion of Zr and Ti atoms to the rotation of the polarization 
in this monoclinic phase. One may also want to know 
if the local (microscopic) structure of a PZT alloy near 
its MPB differs strongly from its average (macroscopic) 



structure, as suggested by Ref. ||. 

Moreoverp-PZT is strongly piezoelectric in the vicinity 
of the MPBd, which is a primary reason why these ma- 
terials are of cucrent use in transducers and other piezo- 
electric devicesQ. Recent theoretical and experimental 
work has reexamined the pieauelectricity of tetragonal 
and rhombohedral PZT alloyailfl. However, we are not 
aware of any previous studies of the piezoelectric or di- 
electric response in the monoclinic phase of PZT. 

The purpose of the present article is to apply the 
newly-developed first-principles derived scheme of Ref. 
^ to PZT and to investigate all the issues mentioned 
above. Our main findings are as follows. (1) As the 
Ti composition decreases in the recently discovered mon- 
oclinic phase, the polarization associated with Zr atoms 
rotates towards the [111] direction and grows in magni- 
tude, while the polarization associated with Ti atoms lags 
in its rotation toward this [111] direction and shrinks in 
magnitude. (2) There are strong fluctuations of the di- 
rection and magnitude of the local polarizations centered 
in the different unit cells of PZT alloys near the MPB, 
indicating that various unit cells adopt a local structure 
somewhat different from the macroscopic average. (3) 
The recently discovered monoclinic phase exhibits very 
large piezoelectric and dielectric coefficients. 

The paper is organized as follows. In Sec. II, we de- 
scribe our approach in detail. Section III reports our 
results, and we conclude in Sec. IV. 
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II. METHODS 

We use the numerical scheme proposed in Ref. ^ which 
consists of constructing an alloy effective Hamiltonian 
from first-principles calculations. For a ferroelectric ma- 
terial, the effective Hamiltonian should include structural 
degrees of freedom corresponding to the ferroelectric lo- 
cal soft mode and the strain variables. These are the 
most important degrees of freedom because ferroelectric 
transitions are accompanied by a softening of the phonon 
soft mode and by the appearance of a strainQ. Moreover, 
a realistic alloy effective Hamiltonian must also include 
the compositional degrees of freedom, because the atomic 
arrangement caji strongly affect the ferroelectric proper- 
ties of an alloylj. We proposed to incorporate all such 
degrees of freedom by writing the total energy E as a 
sum of two energies, 

E{{U^},{V,},^]H{(TJ}) = 
£^avo({uJ, {Vi},r]H) 

+ Sioc({u,},{v,},{a,}) , (1) 

where is the local B-centered soft mode in unit cell i of 
the ABO3 perovskitc material under study; {v^} arc the 
dimensionless local displacements which are related to 
the inhomogeneous strain variables inside each cell and 
which are centered on the A-siteaj; rjH is the homoge- 
neous strain tensor; and the {ctj } characterize the atomic 
configuration of the alloy. That is, aj=+l or —1 corre- 
sponds to the presence of a Zr or Ti atom, respectively, at 
the B-sublattice site j of the Pb(Zri_2;Tia;)03 solid solu- 
tion. The energy Save depends only on the soft-mode and 
strain variables. The {cFj} parameters are incorporated 
into the second energy term Sioc, which thus accounts 
for the chemical differences between Zr and Ti atoms. 

Expresainns for the total energy E have recently been 
proposedBllj for simple ABO3 perovskite systems (i.e., 
in the absence of {(Jj} variables). These have been 
very successful both for reproducing phase transition 
sequettcesEnld'El and for studying ferroelectric domain 
wallalj, as well as for calculating finite-tepu^rature di- 
electric and electromechanical propertiesE^tJ. Here, for 
i^ave, we generalized the analytical expression of Ref. || to 
the case of the Pb(Zri_a;Tia;)03 alloy, bjj-making use of 
the virtual crystal approximation (VCA)ll3il3. We thus 
replaced the Pb(Zri_j;Tij;)03 alloy by a virtual (uniform) 
Pb(B)03 simple system in which the (B) atom is a virtual 
atom involviutt a kind of potential average between Zr 
and Ti atomsO. iSavc thus consists of five parts: a local- 
mode self-energy, a long-range dipole-dipole interaction, 
a short-range interaction between soft modes, an elas- 
tic energy, and an interaction between the local modes 
and local strainil. The analytical expression for i?avc 
has 18 free parameters that are determined by fitting 
to the results of almost 40 first-principles calculations on 
small VGA cells (typically between 5 and 20 atoms/cell) 
following the procedure of Ref. ^. More precisely, the 



first-principles method used is the plane-wave ultrasoft- 
pseudopotential metliodE£l within the local-density ap- 
proximation (LDA)|2l, and the VGA approach adopted 
is the one of Ref. Table I reports the resulting 18 
expansion parameters of -Bavo for the Pb(Zro.5Tio.5)03 
alloy. 



TABLE I. Expansion parameters (in atomic units) of 
for Pb(Zro.5Tio.5)03. 



Parameter 


Value 




0.0138 


a 


0.011 


7 


0.002 


ji 


-0.00577 


ji 


n ni /I op; 

U.U14Z0 


33 


0.00140 


ji 


-0.00094 


js 


0.00141 




0.00006 


h 


0.00003 




5.22 


B12 


1.67 


B44 


1.22 


B±xx 


-0.374 


Blyy 


-0.155 


B^yz 


-0.068 


z* 


7.342 


^00 


7.150 



While expressions were available for iSave , we were not 
aware of any analytical expression that had previously 
been proposed and tested for i^ioc- Following the spirit 
of the "computational alchemy" method developed for 
calcuLalips; the compositional energy of semiconductor 
alloysEJa, we derived Eioc by treating the alloy config- 
uration {aj} as a perturbation of the VGA system. We 
adopted an expression that includes: (i) the on-site effect 
of alloying on the self-energy up to the fourth order in 
the local mode amplitude ; and (ii) the intersite contri- 
butions involving the first-order terms in a perturbation 
expansion in powers of Uj (i.e., terms that are linear in 
Uj or Vj). That is, 

-Bioc({Uj},{v,},{CTj}) = 

i 

+ '^[Q\j-i\ ^3 Oj-j • + R\j-i\ fji • Vj] , (2) 

where the sum over i runs over all the unit cells, while 
the sum over j runs over the mixed sublattice sites. Ui^^ 
Uiy and are the Gartesian coordinates of the local- 
mode Ui. Bji is a unit vector joining the site j to the 
center of the soft mode u^, and fji is a unit vector join- 
ing the site j to the origin of v^. Aa((Ti) and A'j(ai) 
characterize the on-site contribution of alloying. Their 
strength and sign reflect how the identity of the atom 
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sitting on the i site affects the local-mode self energy of 
the "VGA" Pb(Zri_;^Ti^)03 soHd solution. Q\J-^\ and 
R\j-i\ are related to intersite interactions between the al- 
loy parameter aj on the site j and the local mode Ui and 
the strain-related Vi at the site i, respectively. Q\j-i\ and 
R\j-i\ only depend on the distance between i and j up 
to the third neighbors shell, while for symmetry reasons, 
the expression for the intersite interactions becomes more 
complex when going beyond the third neighbors shell. In 
principle, terms involving higher powers of {ffj}, and 
Vi might be included to improve the quality of the ex- 
pansion, but as shown in Ref. ^, we found this level of 
expansion to give a very good account of experimental 
findings. We also found that Q\j-i\ and R\j-i\ rapidly 
decrease as the distance between i and j increases. As a 
result, we included the contribution up to the third neigh- 
bors for Q\j-i\ (denoted Q ^ , Q2 and Qg in the following) 
and up to the first neighbor shell for R\j-i\ (denoted Ri). 

The parameters Aa(CTi), A'y{ai), Q\j-i\ and R\j-i\ 
are also derived by performing first-principles calcula- 
tions, and are given in Table II in the case of the 
Pb(Zro.5Tio.5)03 alloy. More precisely, Aa{(Ji) and 
A'y{ai) are derived by computing the energy of a 5- 
atom cell containing a true B-atom [e.g., Ti or Zr in 
Pb(Zr,Ti)03] when the atoms are displaced as in the 
VGA local mode. Then Q\j-i\ and R\j-i\ are derived 
by using large ideal cubic cells (up to 40 atoms) contain- 
ing a central true B-atom surrounded by VGA atoms, 
and are simply related to the atomic forces occurring on 
these VGA atoms. 

TABLE 11. Expansion parameters of Eioc (in atomic units) 
for Pb(Zro.5Tio.5)03. a—+l and —1 corresponds to the pres- 
ence of a Zr and Ti atom, respectively. 



a + Aa{-1) = 
7 + A7(+l) = 
7 + A7(-l) = 



0.008 + 0.02 (0.5 -x) 
- 0.008 - 0.14 (0.5 -x) 
0.005 - 0.02 (0.5 -x) 



Parameter 


Value 


Aa{+1) 


0.003 


Aa(-l) 


-0.003 


A7(+l) 


-0.010 


A7(-l) 


0.003 


Qi 


0.00160 


Q2 


-0.00028 


Qa 


-0.00018 


Ri 


-0.0125 



(3) 



Such a linear composition-dependence approach is only 
realistic when exploring a narrow range of compositions, 
as done in Ref. | and in the present study. 

Once our effective Hamiltonian is fully specified, the 
total energy of Eq.(l) is used in Monte-Gar lo simulations 
to compute finite-temperature properties of PZT alloys. 
We use 10 X 10 X 10 supercells (5000 atoms), since this 
choice yields well- converged results at low temperatureEJ. 
The {aj } variables are chosen randomly in order to mimic 
maximal compiQaitional disorder, consistent with exper- 
imental realityE3, and are kept fixed during the Monte- 
Garlo simulations. We find that averaging our results 
over a couple of different realizations of the disorder leads 
to well-converged statistical properties. The outputs of 
the Monte-Garlo procedure are the local mode vectors 
u and the homogeneous strain tensor rjH- We use the 
correlation-function approach of Refs [l^jlj to derive the 
piezoelectric and dielectric response from these Monte- 
carlo simulations. Up to 10^ Monte-Garlo sweeps are 
first performed to equilibrate the system, and then 2x 10'' 
sweeps are used to get the various statistical averages. In 
the present study, the temperature is kept fixed at 50 K. 
Note that Refs. demonstrate that our approach leads 
to a (converged) theoretical Gurie temperature Tctheo of 
1032K for Pb(Zro.5Tio.5)03, which is |iHuch higher than 
the experimental value T^^exp of 640 KcS. This difficulty 
of reproducing Tc seems to ber-a^-ftHieral feature of the 
effective-Hamiltonian approachSM^til, and may be due 
to higher perturbative terms not included in the analyt- 
ical expression for the total energy. When comparing 
with measurements, this shortcoming can be overcome 
by multiplying the temperature used_in the simulation 
by a constant factor of Tc^exp/Tcthei^^- As a result, our 
simulated temperature of 50 K corresponds to an "exper- 
imental" temperature around 30 K. 



III. RESULTS 



In principle, all the parameters in Eqs. (1) and (2) 
should change when one varies the composition x in the 
Pb(Zri_a;Tia;)03 solid solution. However, we assumed 
that only the parameters related to the local-mode self- 
energy - i.e., a and 7 for the VGA alloy, and Aa{ai) 
and Aj{ai) in Eq (2) - can significantly change with 
composition. This composition-dependence was assumed 
to be linear, and was determined by performing first- 
principles simulations on cells with two different com- 
positions, namely a;=0.5 and x = 0.45. The resulting 
composition-dependencies are 

a + Aa(+1) = 0.014 + 0.02 (0.5 - x) 



A. Structural properties 

The averaged homogeneous strain variables obtained 
in Pb(Zri_a;Tia;)03 from our simulations are shown in 
Fig. 1 as a function of the composition x. These strain 
variables are measured relative to the LDA-calculated 
minimum-energy cubic structure with lattice constant 
flo = 7.56 a.u., and are expressed in the Voigt nota- 
tion. For Ti compositions larger than 49%, we have 
Vi ^ 12 7^ 0, ?73 > r]2, and 1J4 ~ 7/5 ~ rje = 0. This strain 
tensor corresponds to a tetragonal phase with 5 atoms per 
unit cell. For compositions lower than 47.5%, we predict 
that the system adopts the so-called "high-temperature" 
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FIG. 1. The average homogeneous strain tensor t/h as a 
function of composition in disordered Pb(Zri_i:Tia;)03 at 
50 K. Strains are measured relative to the theoretical mini- 
mum-energy cubic structure of lattice constant 7.56 a.u. 

(S-atoms per unit cell) rhombohedral phase0 since 771 = 
^72 = 773 7^ and 774 = 775 = 775 7^ 0. The most interest- 
ing feature of Fig. 1 is the behavior of the homogeneous 
strain for the compositional range between 47.5% and 
49%: as x decreases, (1) 771 and 772 slightly increase and 
remain equal to each other, (2) 773 strongly decreases but 
is still larger than 771 , (3) 774 and 775 increase and are equal 
to each other, and (4) 77g increases and is smaller than 774. 
This behavior is characteristic of an intermediate phase 
that is neither tetragonal nor rhombohedral, but rather 
is monoclinic. From the strain tensor shown in Fig. 1, we 
further predict that this monoclinic phase can be charac- 
terized by a 10-atom conventional cell (5-atom primitive 
cell) with lattice vectors = ao(— 1 — 771 — 775/2, —1 — 
771-775/2,-774), h„i = 00(1-^771-776/2, -1-771+775/2,0), 
and Cm = ao{r]4/2, 774/2, 1 -|- 773). These predictions are 
in excellent quantitative agreement with the lattice vec- 
tors of the monodinic phase experimentally discovered 
by Noheda et aln. Fig 1 also reveals that the strain 
variables continuously change with composition when 
crossing the tetragonal-to-monoclinic transition, while 
these variables exhibit a sudden jump at the monoclinic- 
to-rhombohedral transition composition. Interestingly, 
these two distinct features are consistent with the predic- 
tions of Ref. ^ i.e., that the tetragonal-to-monoclinic 
transition is of second order while the monoclinic-to- 
rhombohedral transition is of first-order. 

Our simulations also agree with measurements for the 
narrowing of the compositional range otthe monoclinic 
phase when increasing the temperatureo. For instance. 
Fig. 1 indicates that the monoclinic phase is predicted 
to occur for 0.475 < x < 0.49 at r=50K, while we 
predict (not shown here) that increasing the simulated 
temperature up to 485K - which corresponds to an ex- 
perimental temperature of 300 K - leads to the existence 
of the monoclinic phase for 0.475 < x < 0.485. As 
observed in Ref. ^ the tetragonal-to-monoclinic tran- 
sition composition thus decreases when increasing tem- 
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Composition x 

FIG. 2. Average cartesian coordinates (ux), {uy) and (uz) 
of the local-mode vector as a function of composition x in 
disordered Pb(Zri_3;Tia;)03 at 50 K, corresponding to an ex- 
perimental T=30K after rescaling (see text). 

perature while the rhombohedral-to-monoclinic transi- 
tion concentration is independent of the temperature. 
Unfortunately, the calculations are not sufficiently pre- 
cise to determine whether the phase diagram of PZT 
exhibits a triple point where the rhombohedral, mono- 
clinic and tetragonal phases meet at a given composition 
and temperature as suggested in Ref. |29|, or whether the 
monoclinic phase instead survives right up to the bound- 
ary with the paraelectric cubic phase. Furthermore, Ref. 

reports that, at r=20K, the monoclinic phase is ob- 
served for 0.46 < a: < 0.51 , i.e. for a larger range that 
the one shown in Fig. 1. This may be explained by our 
empirical finding that the range of the monoclinic phase 
depends strongly on the parameters of Eq. (3). A slight 
adjustment of these parameters may thus lead to a better 
agreement with experiment for the compositional range 
of the monoclinic phase. Note also that our model can- 
not prcdicvt-|the so-called "low-temperature" rhombohe- 
dral phasaHIl since this involves oxygen octahedra-tilting 
degrees of freedom that are not including in the present 
model. Consequently, our simulations demonstrate that 
the oxygen tilts are not the driving force for the transi- 
tion or for the occurrence of the monoclinic phase. 

Figure 2 shows the cartesian coordinates {{ux), (uy) 
and (uz)) of the supercell average of the local mode vec- 
tors in Pb(Zri_a;Ti2,)03 as a function of the composition 
X atT — 50 ii', as predicted by our approach described by 
Eqs. (1) and (2). The average local mode (u), and hence 
the polarization, is parallel to the pseudo-cubic [001] di- 
rection in the tetragonal phase {x > 0.49%) while the 
polarization becomes parallel to the pseudo-cubic [111] 
direction in the rhombohedral phase {x < 0.475%). In 
the monoclinic phase, the change of strain shown in Fig. 1 
is associated with a decrease of (uz) while (ux) and (uy) 
increase and remain nearly equal to each other as the 
composition x decreases from 49% to 47.5%. Figure 2 
thus demonstrates that the electrical polarization rotates 
from the pseudo-cubic [001] direction to the the pseudo- 
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FIG. 3. Angle between the average local-mode vector and 
the pseudo-cubic [001] direction as a function of composition 
in disordered Pb(Zri_a;Tij;)03 at 50 K. Circles and diamonds 
refer to averages over local modes centered on Zr and Ti sites, 
respectively; solid line refers to average over all local modes. 
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FIG. 4. Magnitude of the average local-mode vector as 
a function of composition in disordered Pb(Zri_a;Ti3;)03 at 
50 K. Circles and diamonds refer to averages over local modes 
centered on Zr and Ti sites, respectively; solid line refers to 
average over all local modes. 



smaller than the magnitude of uxi in the entire mono- 
clinic phase concentration range, and that the magnitude 
of uzr increases while the magnitude of uxi decreases 
when the Ti composition decreases from 49% to 47.5%. 
This leads to a nearly composition-independent magni- 
tude of (u). In other words, the total polarization simply 
rotates in the monoclinic phase while the polarization as- 
sociated with Zr (respectively, Ti) atoms rotates and also 
grows (respectively, shrinks), as the Ti concentration de- 
creases. 

The authors of Ref. || suggested that near the MPB, 
the rhombohedral and tetragonal phases of PZT can be 
described in terms of a structure that is locally mono- 
clinic, and in which the average polarization along the 
[111] or [001] pseudo-cubic direction can occur by means 
of fluctuations between a subset of three or four nearby 
monoclinic orientations. In this picture, the transition 
from the rhombohedral or tetragonal phase to the mon- 
oclinic phase would occur by the freezing in of one of 
these monoclinic orientations. Our effective Hamiltonian 
can be used to investigate the local structure of PZT 
around its MPB, and, in particular, to check if the local 
structure of the tetragonal or rhombohedral phase is dif- 
ferent from its average structure. Figures 5(a-c) display 
the predicted local modes distributions, at T=50K, for 
Pb(Zri_a;Tia;)03 solid solutions with a;=0.50 (tetragonal 
average structure), a;=:0.482 (monoclinic average struc- 
ture) and a;=0.47 (rhombohedral average structure), re- 
spectively. One can clearly see that the direction and 
magnitude of the local modes fluctuate around their av- 
erage value in any of these three PZT solid solutions. As 
a result, various unit cells adopt a local structure different 
from the macroscopic average. However, our calculations 
do not support the hypothesis that the tetragonal and 
rhombohedral phases of PZT are simply made up of a 
small number of local monoclinic phases, since Figs. 5(a) 
and 5(c) do not show the distributions breaking up into 
clusters centered along the monoclinic directions. 



cubic [111] direction as the Ti composition decreases in 
this monoclinic phase. 

To better understand the separate contribution of Zr 
and Ti atoms to this rotation. Fig. 3 shows the composi- 
tion dependence of the angle S between the pseudo-cubic 
[001] direction and the uzr (respectively, uxi) local mode 
vectors averaged over cells centered on Zr (respectively, 
Ti) atoms. The angle between the pseudo-cubic [001] di- 
rection and the entire supercell average (u) is also shown 
in Fig. 3. Similarly, Fig. 4 displays the variation of the 
magnitude of (u), uzr and uxi as a fimction of the com- 
position. Figure 3 demonstrates that, in the monoclinic 
phase, uzr and uxi also rotate from the [001] direction 
- for which S = - to the [111] direction - for which 
S = 54.7° - as the Ti composition decreases from 49% 
to 47.5%. In this monoclinic phase, the pseudo-cubic 
[111] direction is always closer to uzi than to uxi. Fur- 
thermore, Fig. 4 indicates that the magnitude of uzi is 



B. Piezoelectric and dielectric properties 

We now use our alloy effective Hamiltonian to in- 
vestigate the piezoelectric and dielectric properties of 
Pb(Zri_:rTia;)03 vs. a; at T = 50K. Figure 6 shows the 
piezoelectric coefficients (^33 and dis as a function of the 
composition a;, when representing the piezoelectric ten- 
sor in the orthonormal basis formed by ai — [100], a.2 = 
[010] and as = [001]. Figure 7 displays the composition 
dependence of the dielectric susceptibilities xii and X33 
expressed in the same basis. 

One can notice that dis is much larger than (^33 in the 
tetragonal phase, i.e., for x > 49%. This is consistent 
with recent measurements revealing that the piezoelec- 
tric elongation of the tetragonal unit cell of PZT does not 
occur along the polar [001] direction!]. The large value 
of di5 also explains the strong piezoelectric response ob- 
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FIG. 5. Projection of the distribution of local-mode orien- 
tations for disordered Pb(Zri_3;Tia;)03 at 50 K for (a) a;=0.50, 
(b) a;=0.482, and (c) x=0A7. 
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FIG. 6. Piezoelectric coefficients dsa and dis as a function 
of composition in disordered Pb(Zri_a;Tii:)03 at 50 K. Statis- 
tical errors are estimated to be ~10% of the values displayed. 
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served in ceramic samples, since this latter involves 
average over the single- crystal coefficients di^ and 1^331 

One can also see that ^33 is predicted be quite large 
in the rhombohedral phase (occurring for x < 47.5%). 
This prediction agrees with the theoretical and experi- 
mental findings of Ref. |3^:^ that the ^33 coefficient of 
rhombohedral materials can be very large along the [001] 
direction, i.e. away from the polar direction which is 
oriented along the pseudo-cubic [111] direction. 

Fig. 6 also reveals that 1^33 reaches its largest value near 
the monoclinic-to-rhombohedral transition, while dis 
peaks near both the monoclinic-to-rhombohedral and 
tetragonal-to-monoclinic transitions. The most striking 
feature of Fig. 6 is that dis has a remarkably large value 
- above 600 pC/N - in the entire monoclinic phase range. 

Fig. 7 demonstrates that the dielectric susceptibility 
X33 behaves in a similar way as the piezoelectric coeffi- 
cient dss, in the sense that X33 is also peaked near the 
monoclinic-to-rhombohedral transition and is also much 
larger in the rhombohedral phase than in the tetragonal 
structure. On the other hand, xii behaves in a manner 
similar to dis, since they both have peaks at both tran- 
sitions, and are both very large in the monoclinic phase. 

It thus appears that the rotation of the polarization 
(see Fig. 2) not only leads to a very high ^15 piezoelectric 
response (see Fig. 6), consistent with the finding of Ref. 
[3^ , but also to a large xii dielectric response (see Fig. 7). 



IV. CONCLUSIONS 

In summary, we have used the first-principles derived 
computational scheme proposed in Ref. ^ to study low- 
temperature properties of disordered Pb(Zri_j,Tij,)03 
solid solutions near the JVIPB. 

We find that the monoclinic phase reported in Ref. ^ 
acts as a structural bridge between the tetragonal phase 
{x > 49% at T=50K) and the rhombohedral phase 
{x < 47.5% at T=50K), with the electric polarization 
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FIG. 7. Dielectric susceptibilities xii and X33 as a function 
of composition in disordered Pb(Zri_a;Tia:)03 at 50 K. Statis- 
tical errors are estimated to be ~10% of the values displayed. 



rotating from the pseudo-cubic [001] to [111] direction as 
the Ti concentration decreases through the monoclinic 
range. The polarization associated with the Zr atoms 
differs from the polarization associated with Ti atoms in 
the monoclinic phase, in that the former rotates faster to- 
wards the [111] direction and grows in magnitude, while 
the latter rotates more slowly and shrinks as the compo- 
sition X decreases. Furthermore, we also investigated the 
local structures of tetragonal, monoclinic and rhombohe- 
dral PZT alloys, and found that the local polarizations 
centered in the different unit cells fluctuate significantly 
in both magnitude and direction around their average 
value. However, the pattern of these variations does not 
support a picture in which the tetragonal or rliombo- 
hedral phases could be regarded as arising from fluctua- 
tions among neighboring monoclinic states. Finally, some 
piezoelectric and dielectric coefficients are predicted to be 
extremely large in this monoclinic phase. 
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